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Excitation of plasmons in two-dimensional electron gas with defects by microwaves: 

Wake-field method 

Eduard Takhtamirovt] and Roderick V. N. Melnik^'^ 

^M^NeT Laboratory, Wilfrid Launer University, Waterloo, Ontario N2L 3C5, Canada 

^BCAM, Bizkaia Technology Park, 48160 Derio, Spain 

We develop an analytical method to find plasmons generated by microwaves in a two-dimensional 
electron gas with defects. The excitations are expressed in terms of the wake field of a charged 
particle moving in plasma. The result explicitly addresses the efficiency of the photon-to-plasmon 
conversion and the type of excitation. While strong absorption of the radiation by the excitations 
is reached at larger plasmon wave numbers, intense persistent plasma waves are created at optimal 
ones. The latter wave numbers depend on the spectrum of plasmons and the distance that the waves 
are required to travel without being substantially attenuated. Their type, which can be traveling 
or standing, is governed by the geometry of the defects and the polarization of the radiation. We 
identify such types of traveling plasmons as circular plasmons, excited at dot defects, and traverse 
plasmons, excited at straight wire defects and traveling away from (toward) the wires if their group 
velocity is positive (negative). Nonlinear excitations are also accounted for. In particular, we analyze 
the zeroth harmonic, linear in the microwave intensity, which has the character of a frozen charge 
density wave. The interference of elementary wakes from defects arranged in truncated periodic sets 
can produce amplified plasmons, which are easily portrayed. 

PACS numbers: 73.20.Mf, 78.70.Gq 



I. INTRODUCTION 

The problem of excitation of plasmons has always re- 
ceived considerable attention associated primarily with 
experimental verification of theoretical predictions for 
spectra of quasi-particles in interacting many-electron 
systems J^^ Recently, the interest has been greatly fueled 
by the development of a new generation of photonic de- 
vices, with plasmons playing a key roleii"— The efficiency 
of the photon-to-plasmon conversion is always a major 
concern. To generate intense plasma waves, one has to 
deal with a conflict between the requirement to have a 
system with high electron mobility, beneficial for plas- 
mons, and the necessity to build in a strong enough po- 
tential that breaks the system's continuous translational 
symmetry, inevitably impairing the mobility. Different 
balanced schemes have to be adopted depending on the 
specifics of the task^^^ - 

The type of excited plasmon is also an important sub- 
ject connected to the destination of the excitations, or 
explicitly to the mechanism of their terminal detectiortii 
or transformation.^ If a particular situation requires the 
excitation of, for example, traveling plasmons, a stand- 
ing charge density wave should be treated as detrimental 
rather than just idle because merely the presence of its 
exciting agent deteriorates the whole system's character- 
istics. This leads to the efficiency issue again. 

For defects arranged in a grating with a spatial pe- 
riod 5x, it was found that light can effectively excite 
plasmons with wave numbers Q — iim/Sx, where n — 
l,2,...i^"'* This gives only a qualitative picture lack- 
ing important details. A deeper insight usually requires 
a direct solution of Maxwell's equations. This can be 
performed (semi) analytically while considering surface 
plasmon-polaritonsi^"— in metals, by using local dielec- 



tric constantS j^^d^ or in two-dimensional (2D) semicon- 
ductor systems still within the approximation of local 
conductivities)^ which does not allow description, in par- 
ticular, of collective excitations with large wave numbers. 
The Green's-function-based approaches are probably the 
only means to deal with such situations as the excitation 
of plasmons by spatially inhomogeneous electromagnetic 
beamai^iii or scattering of the plasmons themselves at 
defects j "^' ' '-'^' But such integral equation methods would 
invoke very cumbersome calculations with an obscure re- 
sult if they were applied to systems^ exhibiting a signifi- 
cant space-time dispersion of the response function. 

We develop an effective alternative approach appli- 
cable in situations where the homogeneous microwave 
(MW) radiation acts on an electron gas with defects, 
which can be treated as weak scatterers. It is analo- 
gous to the formalism of the polarization wake of a fast 
charged particle in plasmai ^^i^^ Here it is the static de- 
fect, which plays the role of a charged particle, that ex- 
cites a wake in the electron gas, the latter moving as a 
whole under the action of the MW field. For simplicity, 
we put two restrictions on the system under study. First, 
we regard collective excitations for which the retarda- 
tion of the electromagnetic interaction can be neglected. 
In particular, we focus our attention on a 2D electron 
gas (2DEG) with plasmons^^ or magnetoplasmons^^iSi 
(in quantizing magnetic fields) as the lowest-lying exci- 
tations. Actually, surface plasmons with large wave num- 
bers, for which the admixture of polaritons is negligibly 
small, are also inferred. We refer to all these excitations 
as plasmons unless otherwise stated. And second, we 
suppose that the internal structure of the defects is not 
changed by the MW field, which affects only the 2DEG. 

The expression we derive below for excited plasmons 
is then analyzed in terms of the efficiency of the photon- 



to-plasmon conversion and the type of the resulting exci- 
tations. Our results indicate that generally two goals of 
the conversion should be distinguished from one another: 
the generation of intense persistent plasma waves and the 
attainment of strong absorption of the radiation by the 
excitations. While the former is crucial for fast plasmonic 
devices^ the latter is commonly used to study the dis- 
persion of plasmons as suchri^— We show that maximum 
efficiencies of these two processes are generally reached at 
different plasmon wave numbers and thus different values 
of the radiation frequency, depending on the spectrum of 
plasmons. 

The paper is organized as follows. In Sec. HI] we re- 
view general properties of a 2DEG, the single-electron 
spectrum of which is parabolic, with embedded defects 
under the action of MW radiation with a homogeneous 
electric field. In Sec. Illll we focus on the derived expres- 
sion for the screened potential of the defects and specify 
the plasmon contribution. Wakes from a wire defect and 
a truncated periodic set of wire defects are dealt with 
in Sees. |IV] and El respectively. In Sec. |VI] we study 
wakes from dot defects. In Sec. lVIII we discuss the effects 
of retardation and the MW-induced polarization of the 
defects, neglected in previous sections. Conclusions are 
given in Sec. IVIIII 



II. ELECTRON GAS WITH DEFECTS IN MW 
FIELD 

We consider an isotropic 2DEG unbounded in the 
(a;, y) plane of Cartesian coordinates. The system is ir- 
radiated by MWs with frequency 51 > and in-plane 
electric field E(f) — {Ex{t),Ey{t)) with the amplitudes 



E,{t)=S,e-''''+S*e'''\ j ^ x,y, 



(1) 



where £* is the complex conjugate of £j. The Hamilto- 
nian of a 2D electron with charge — e < and effective 
mass m* is 



Ho 



2m 



-+eE(t)r, 



(2) 



where p — (px,Py) is the momentum operator, which 
includes the vector potential of the external magnetic 
field (if present), and r = {x,y) is a 2D coordinate. 

It is convenient to change the reference frame to the 
one associated with the instantaneous position of the 
classic electron or, in magnetic fields, of the classic cen- 
ter of the electron's cyclotron orbit. For the Schrodinger 
equation, such a transformation, realized with a unitary 
operator S, was proposed more than five decades ago.^^ 
Recently, it has been used by a number of researchers 
(see the work by Dmitriev et alM- and references therein) 
in connection with the experimentally observed MW- 
induced resistance oscillations (see the pioneering work 
by Zudov et al.)^ 



Below, we use two properties of this transformation. 
First, the transformed Hamiltonian Hq does not contain 
the radiation field in this moving reference frame and co- 
incides with the Hamiltonian for the rest reference frame 
without radiation: 



Hq — S Hq S 



2m* 



(3) 



where S^ is the Hermitian conjugate of the operator S. 
And second, the transformation is homogeneous, the co- 
ordinate being transformed as follows: 



?:=S'trS' = r-ro(i), 
where ro(t) = {xo{t),yo{t)) with 

yoit) =re-''"*-f r*e*"*, 



and the amplitudes are 



X = 



e {^£x — i^c£y) 
m*n{ujl-9?) '' 



Y = 



e {Vl£y + iujc£x) 
m*fl (lu^ 



n^ 



(4) 



(5) 



(6) 



Here ujc is the cyclotron frequency, which is zero in the 
absence of an external magnetic field. 

Let us consider now an ensemble of electrons in thermal 
contact with a reservoir. If we apply the above coordinate 
transformation to each electron and each microscopic el- 
ement of the reservoir, the many-electron Hamiltonian 
will also lose the information about the external radia- 
tion. This is so because the homogeneous transforma- 
tion does not change the inter-electron pair interaction. 
On the other hand, physical properties of the reservoir, 
which has infinite degrees of freedom, will not be affected. 
Therefore, in the moving reference frame we can intro- 
duce an electron response function coinciding with the 
one for the same system in the rest reference frame with- 
out radiation. 

As the next step, let the 2DEG have a single defect, 
treated as a weak perturbation, with the 2D potential 
Vb(r) in the rest reference frame. We suppose that its 
internal structure is not developed under the MW field, 
the defect thus resembling an elementary particle. In 
the moving reference frame it is only its position that 
depends on time: 

Fo(r, t) = S^Vo{r)S = V„{r - r„{t)). (7) 

For the 2D space Fourier transform of the above we have 

Fo(q,t) = Fo(q)e-"»"-«(*\ (8) 

where Vq (q) is the Fourier transform of the potential of 
the defect in the rest reference frame and q = (qx,qy) is 
a 2D reciprocal coordinate. The time Fourier transform 
of Eq. (O yields 

-t-oo 

%{ci,uj)^Vo{ci) J2 A,,{q)S{nn-oj), (9) 



where ui is the frequency domain variable, 5{uj) is the 2DEGs, similar to the case of hot electrons in an electron 
Dirac delta function, and for A„ (q) = An{qx,qy) we gas.^® 



have 



27r J 



Restoring the direct space and time coordinates for 
Eq. pT|) . we have 



(10) 



-ir/n 



V{v,t)^ J2 d\ 



+°° " ,2 yo(q)A.(q)e^'»""-"''* 



so that A* (q) = yl_„ (— q)- The potential of the defect 
screened by the 2DEG is 



e{q,nn) 



(12) 



F(q,c^) = 



^o(q) 
e{q,uj) 



+ 00 



Y^ An{q)5{nVL- 



(11) 



n— — oo 



In the rest reference frame the screened potential of the 
defect is 



y(r, t) = SV{v, t)S^ = V{v + ro(t), t). (13) 



where e (g, u) is the permittivity of the 2DEG without 
radiation. Here, we do not account for exchange interac- 
tion processes, which exist when the defect contains lo- Reduced by its bare value, the potential is 5V{Y^t) 
calized electrons. They may be important for low-density V{y, t) — Vb(r, t): 

I 



5V{v, t)= Y. j ^'9 ^o(q) i^Ul^) - l) ^" (^) ^"' ("^) e^'i--'("+"')"*. 



(14) 



n, n — — oo 



The remaining part of this paper is devoted to the anal- 
ysis of Eq. ([T4l). We consider two elementary types of 
defect: a straight wire breaking the translation symme- 
try in the x direction or a dot breaking the translation 
symmetry in both in-plane directions. 



III. EXTRACTION OF PLASMON FIELD 



time r > 0: 



Equation ([14]) resembles the formula for the screened 
potential of a charged particle moving in plasma.— "^i In 
particular, the radiation harmonics nVi, n = ±1,±2, . . ., 
play the role of the Doppler frequency shift Slj ~ qv 
in the argument of the response function, where v is 
the velocity of the charged particle. The roots of the 
equation e (q, nfl) = define the wake field or, in other 
words, plasma excitations, characterized by the wave vec- 
tor Req — Q. The harmonic n = in Eq. p^ does 
not relate to the excitation of real plasmons merely con- 
tributing to the "static" screening, not considered in this 
paper. 

We suppose that the function Vo(q) has no poles on 
the real axis, the oscillatory behavior of the screened po- 
tential thus entirely originating due to the plasma excita- 
tions. This allows us to calculate the oscillatory part of 
the potential analytically without numerical evaluation 
of the integral in Eq. ([U]), using the plasmon pole ex- 
pression for the function e~^ (qj^) near the plasma res- 
onances, with damping characterized by the relaxation 



e{q,uj) 



U) 



uj — ujp{q) + i/r uj + ujp{q) + i/r 



(15) 

Here ujp{q) > represents the spectrum of plasmons, and 
we suppose that the electron gas is clean enough so that 

nT» 1. 

Actually, the permittivities of the electron gas in the 
moving and the rest reference frames coincide only in the 
limit 1/t — >■ -|-0 in the absence of disorder. We have sup- 
posed that Eq. p^ still holds for finite 1/t if in the rest 
frame the reciprocal permittivity has a similar form but 
where damping is characterized by the relaxation time 
To > 0. If the approach is valid, generally t 7^ tq. We can, 
however, establish conditions in which t ~ tq to guaran- 
tee the applicability of Eq. P3|) . For definiteness, we re- 
gard disorder due to impurities with the potential W{r), 
which transforms into the potential W{r,t) = S''W{r)S 
in the moving reference frame. For the space-time Fourier 
transform of the latter, analogously to Eq. ([9]), we have 



+00 



W{q,u) = W{q) J2 An{q)5inn 



(16) 



n— — 00 



where iy(q) is the Fourier transform of W{r). In 
Eq. (ITB)) . the harmonics with n ^ define an effective 
"additional" time-dependent scattering potential, whose 
contribution to the relaxation time can be taken into 
account by using the technique developed for electron- 
phonon scatteringJS Nonetheless, we can estimate its 



contribution here and determine when it can be ne- 
glected. 

Let us intr oduce the characteristic amphtude R = 
A/|yp + |Xp, where X and Y are given by Eq. ©. For 
wave numbers q <^ 1/R we have 

Ao{q) = l + 0{q'R'), A„(q)=0((gi?)H). (17) 

As a resuh, due to the scalar character of the scattering 
cross section, which does not depend on the phase of the 
MW field E(t), we can assume that 



2kp, and 



r = To(l + 0(<7f„,pi?2)), 

where for impurities inside the 2DEG gim 
A;f is the electron Fermi wave number, so that we require 
2k-pR ^1. If the scattering is governed by the pres- 
ence of charged impurities separated from the 2DEG by 
a spacer of width Limp > 1/2A;f, we have gimp ^ 1/Limp, 
and the approximation r sa tq is reached if i? ^ Li^p. 

To estimate the above quantities, as an example, we 
consider a GaAs/AlGaAs-based 2DEG with concentra- 
tion A^s = 3 X 10"'^^ cm~^, radiation with angular fre- 
quency SI — 10^"^ s~^, and the electric field £ — 2>Q V/cm 
(^ 1 W/cm^ of radiation energy flux) . We have R ^ 1 K 
while l/ZcF ~ 70 A and typical widths of the spacers are 
hundreds of angstroms,'^'' which results in the estimate 
qf R^ < 10"'^. As another example, we can mention 
that typical values of the MW radiation frequency and 
the strength of the field used in experiments on MW- 
induced resistance oscillations are O ~ (lO-'^-'^-lO"'^^) s~^ 
and f ~ 1 V/cm (see, e.g., Ref. [3l|), so that again 
i? ^ 1 A. These cases demonstrate that the approxi- 
mation T sa To can indeed be valid for rather realistic pa- 
rameters. This guarantees the applicability of Eq. ([T5|) . 
which is still questionable otherwise \i R^ l/2kp. 

We therefore consider QR ^ 1. For example, in a 
GaAs/AlGaAs-based 2DEG with concentration Ng — 



3 X 10 cm ^, frequency fi — 10 



and field 



£ = 30 V/cm, the MW field excites 2D plasmons with 
Q « 2 X 10^ cm-i (see below), so that QR - 10^^. Due 
to Eq. (fT7|) . the principal oscillatory part of SV{r,t) [see 
Eq. P^ ] is enclosed in the potential 



1/W(r,i)= J2 /^o(q)f^^-l 



(19) 



An (q) e^'5'-™"* d^q 



where we have used Ao{q) « 1. From Eq. (fT4|) we can 
also extract a frozen wave potential, smaller than the 
principal contribution V^^'{r,t) by the parameter QR: 






(20) 



|v4„(q)|Vi"-d25. 



Other harmonics of the excited plasmons can be sepa- 
rated analogously. 



IV. WAKE FROM WIRE DEFECT 

Let us obtain the excited plasmon field from a straight 
wire defect placed at a; = 0, for simplicity having a sym- 
metric space potential, so that its Fourier transform is 
Vb(q) = 6{qy)Viu{qx) with V^{-qx) = Ku(<7x)- 

The procedure of evaluation of the right-hand sides of 
Eqs. dm) and ([20]) can be illustrated by considering this 
model integral, corresponding to the terms with n = 1 in 
Eqs. dHl) and dini): 



+ 00 



(18) I{x)^ / dq^V^iq^) 



1 



e(fe,0,f]) 



-1 Mi(g„0)e" 



(21) 
where e {qx, 0,il) — e {q, il) \q^=o- We will estimate it for 
the far zone \x\ > a + 2|A|, where 2a > is the diameter 
of the wire (see Appendix \^ . 

For positive x > a-|-2|A|, the original integration con- 
tour may be closed at Im q^ — > +oo by the half-circle path 
l^a;! = const — >■ CXI, Imqx > 0. The integral over this path 
is zero^ because the expression in parentheses in Eq. ((2T|) 
goes to zero for complex q^ -^ oo. Thus, the right-hand 
side of Eq. (pij) is entirely defined by the poles of the 
integrand, the oscillatory behavior originating due to the 
poles of the function 1/e {qx,0, ^)- To find them, we can 
use Eq. (fT5|) . These poles coincide with those of the func- 
tion l/{fl-u)p{qx, 0)+i/T), where LOp{qx,0) = Wp(g)|,^=o, 
because ft + ujp{qx,0) ^ due to our adopted choice 
fi > and ujp{q) > 0. And vice versa, calculating the 
terms with n = — 1 in Eqs. P^ and (PO)) . we deal with 
the poles of the function l/(—n+u)p{qx, 0)-|-i/r) skipping 
the function l/(— O — ijjp{qx, 0) -I- ilr). 

Continuing with Eq. (pij) . we suppose that the poles 
qx — I'j + iQj (j — Ij 2, . . .) of the integrand are simple, 
with small imaginary parts: \q';'\ <^ \q'A, where the real 
parts g' are solutions to the equation Ljp{qj, 0) = fl, while 



1': = 



1 



^'piq-,o)r' 



(22) 



with uj' {qx, 0) = diijp{qx, 0)/dqx- A general case is briefly 
discussed in Appendix |B] The poles are inside the con- 
tour of integration if a;p(g^, 0) > 0. Accordingly, keeping 
only the oscillatory contributions from the poles, we have 
for positive x > a + 2\X\: 



Tix) 



^'piq'pO) 



E 

J ("J,>0) 



^iqjX-q- x 



(23) 



Note that I [x) — >■ as a; — )• -l-oo due to the exponential 
decrements exp (^—q"x) with q'' > [see Eq. (P^ j. 

Similarly, for negative x < — a — 2|A|, with the original 
integration contour closed at Img^; — )■ — oo, only terms 
with Wp(g^-,0) < should be retained: 



X (x) ^ infl 



E 

jK<o) 



VUq'j)Ai{q'j,0) 

^'piqpO) 



^tq^x-q. X 



(24) 
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FIG. 1. Spectrum of 2D plasmons as a function of qx at 
qy = Q (full line) and the radiation harmonic uj = Q. (dashed 
line) showing two points of intersection at the wave numbers 
q'l and ija = -gi- 



Again, I{x) -^ Q as x -^ -oo for q" < [sec Eq. ([^ ]. 

We limit ourselves with the condition of only two 
plasmon wave numbers contributing to the wake: q'^ = 
~92 = *3 > 0. This case is realized for 2D plasmons 
(see Fig. [T]), having the spectrum uip{q) = ui2d{q), with 
^2D{q) = \/27r/Vse^g/m*K, where n is the permittivity of 
the surrounding medium, generally a function of q.^^ It 
holds also for magnetoplasmons, at wave numbers much 
smaller than the reciprocal Larmor radius, with the spec- 
trum ujp{q) ~ \J^1 + uj'^ijiq). The same can take place 
for surface plasmons, depending on the material and the 
spectral rangeJ^ We have plasmons with positive group 
velocity: 



c^;(q'i,o) = -c^;(g^,o) = c^;(Q,o)>o. 



(25) 



For the principal plasmons [see Eq. ([T9l) ]. wc obtain the 
following expression, valid in the far zone |a::| > a + 2|X|: 

,.(1)/ ,x VwSgn jx) cos (Q|a;| - fit + t/j) 

^^'^^ exp(|x|/c.;(Q,0)r) ' ^''^ 

where the undamped amplitude of the potential of the 
plasmons is 



^ w — 



2TTV^{Q)nQ\X\ 



(27) 



and we have used the approximation Ai (g^;, 0) ~ —iqxX, 
while ^ = argX. The phase velocity of these plasmons is 



Vp {x, t) = 



sgn {x) n 



Q + cot (g|a;| - Qt + tP) /lj'JQ, 0)i 



(28) 



They are traveling in the x direction, away from the wire 
on average [neglecting damping, Vp{x) = agn{x)fl/Q]. 
The case when plasmons with negative group velocity 
are excited, resulting in plasma waves traveling toward 
the wire defect, is given in Appendix IBJ 
The frozen wave, derived from Eq. (j20[) . is 



V^'^ix) 



7^^Q|X|sinQ|x| 
exp(|a;|/w;(Q,0)T) 



(29) 



valid in the far zone |a;| > a + 4|X|; see Appendix [K\ Its 
amplitude has an extra small parameter Q\X\ as a factor. 



To determine the structure of the potential P^, , let us 
assume that the wire defect is small on the scales of the 
plasmon wavelength, which means that 2aQ <C 1. In this 
case we can neglect the dependence on Q in the Fourier 
transform Ko(Q): 



a 

VUQ) -^ I v^[x)dx = -(K,), 

27r J TT 



(30) 



where we have introduced the average potential {Vw) of 
the wire defect. Considering, for example, 2D plasmons 
in a medium with a constant permittivity k, so that 
Wp((7, 0) — ujp{q, 0)/2g, we have 



Vr, = 4{V^)a\X\Q^ 



(31) 



Without magnetic fields, X ex l/ft^ [see Eq. ^], so 
Vw oc Q for a fixed radiation intensity. This means that 
without damping it would be more favorable to have ex- 
cited plasmons with larger wave numbers. Exceptions 
exist, obviously, for potentials Vw{Q), not small on the 
scales of the plasmon wavelength, decaying fast enough 
with Q. 

The damped amplitude for 2D plasmons is 



Vn 



■p^^-\MHiQ..O)r ^ Q^-2\MQ/^AQ,0)r^ (33) 



Being a function of the wave number, it reaches a max- 
imum at Q/ujp{Q,0) = T/|a;|. Similar results can be ob- 
tained also for 2D magnetoplasmons and surface plas- 
mons, with specific values depending on details of the 
spectra of the excitations. 

Note that the photon-to-plasmon conversion is still 
more efficient for plasmons with larger wave numbers, 
even though they might not be well sustained in the sys- 
tem due to damping. Thus, corresponding to the stop- 
ping power in the physics of interaction of fast particles 
with plasma , ^"i^"'^ the dissipative current in a 2D electron 
system with impurities in quantizing magnetic fields and 
weak dc electric fields is governed by the excitation of 
plasmons with the largest wave numbers of the order 
of 2A:f'^ The concept of the efficiency of the photon- 
to-plasmon conversion accordingly depends on the goal 
set: to excite intense persistent plasma waves, or to at- 
tain strong absorption of the MW radiation by plasmons, 
which are excited at the defects. Generally, these goals 
are reached at different plasmon wave numbers and thus 
different values of the MW frequency. 

There is a plain consequence of the above, important 
for the problem of the detection of the plasma waves. 
If a plasmon detector^^ is located at a distance L from 
a source of plasmons, it is the most efficient for detec- 
tion of excitations having such wave numbers Q that 
l/ujp{Q,0) ~ t/L. For example, for 2D plasmons this 
condition is Q/ujp{Q, 0) = t/L. Plasmons with much 
larger values of l/ijjp{Q, 0), which normally correspond to 
larger Q, can not arrive at the detector due to damping, 
while plasmons with smaller wave numbers are poorly 
excited. 



Finally, let us estimate the distance L for a 
GaAs/AlGaAs 2D electron system, assuming r = 5 ps, 
which was measured in a system with moderate electron 
mobility,— and using the same parameters as in Sec. IIIII 
Vt = ujpiQ,0) = 10" s-i and = 2 X 10^ cm"!. We 
have L = 2.5 /im or LQ/2t: « 8 plasmon wavelengths. 



V. WAKE FROM TRUNCATED PERIODIC SET 
OF WIRE DEFECTS 



For simplicity, in this section we neglect the finite size 
of the near zone formally considering the limit Q{a + 
4|A|)^0. 

Let us arrange N identical wire defects in a trun- 
cated periodic set (grating) with a period dx = A2t:/Q, 
which is specified by the dimensionless parameter A. 
The values A = 1,2,... provide the plasmon resonance 
condition^"— The functions 

y^") (x, t)^J2 ^^"'' (^ - 27rnA/g, t) (33) 



describe the net potentials of the principal plasma wave 
{m — 1) and the frozen plasmon (m = 0, for which the 
time argument is superfluous). 

The sum in Eq. ([55]) can easily be calculated outside 
the grating region, for x ^ [0,27r(A^— 1)A/Q]. Neglect- 
ing damping for simplicity, we obtain that the amplitude 
of the net plasmon field is proportional to the parameter 
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FIG. 3. Integration contour for Eq. (|39p over complex q for 
the case cos (cf) — 8) > X > 0. The pole P of the function 
e~^ (g, Q) inside the contour defines a plasma wave. 



At A = 1,2,..., we have a = N, resulting in A''-times 
amplified plasmons outside the grating region: 

VJ;}'' {x, t) - -NV^sgn {x) cos {Q\x\ - nt + t/j) , (35) 



V}°\x) ^ NV^Q\X\ sin Q\x\ 



(36) 



(34) 



For non-integer A, we have: \a\ < N. The destructive 
interference cases, for which cr = 0, are obtained for non- 
integer A = n/N with an integer n y^ 0. 

To demonstrate the above, in Fig. [5] we plot the po- 
tential Vj^f (x) from a grating of A^ = 10 wire defects. 
Presented are the constructive interference case A = 1 
and several other arrangements including two destruc- 
tive interference cases at A = 1/2 and A = 9/10. 



VI. WAKES FROM DOT DEFECTS 
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FIG. 2. Potential (in units of VwQlXl) of the net frozen 
plasmon from A = 10 wires, which form a truncated periodic 
grating occupying the space < x < {N — 1)A (in units of 
2tt/Q), for (a) A = 1, (b) A = 3/4, (c) A = 1/2, and (d) 
A = 9/10. Damping is neglected. 



Established for the wire defects, the conclusions on the 
efficiency of the photon-to-plasmon conversion are valid 
for a defect of an arbitrary shape. It is the different 
plasmon wave pattern from a dot defect that makes us 
pay special attention to this case. In addition, the wake 
potentials obtained below can play the role of the Green's 
functions to deal with inclusions of complex shapes. 

Let us focus on a dot defect, placed at x = 0, y — 
0, that has an isotropic potential: Vb(q) = Vo{q). To 
calculate the integrals in Eqs. (|T9)) and ((20)) . we analyze 
the following model integral: 

r (r) = I d'q Vo{q) (^^^ - l) A, (q) e^'i"". (37) 

Using polar coordinates, we transform both the space 
variables x — r cos 6, y — r sin 6, and the integration 
variables q^ — qcoscj), qy ~ qsincf). Being transformed, 
Eq. dniD reads: 






2n 




(38) 



where $ = $ I 



+ CXD 



* = qVo {q)Ai {q COS (/), g sin (p) 

—}—- - \\ e^i^^°<'t'-o) dq. 
£{q,n) J 



(39) 



We will estimate T (r, 6') for the far zone r ^ Ai (see 
Appendix O. In Eq. ^, for -7r/2 + 6* + A < < 
7r/2 + — A, the original integration path pi, which is 
< Reg < i^ with Imq = and real K — > +oo (see 
Fig. 13]) may be supplemented by the quarter-circle path 
P2, which is \q\ = K and < argg < 7r/2, and by the 
path p3 over the imaginary axis < Im q < K to form a 
contour. 

The integral over the path p2 is zero^^ if 



1 



e{q,uj) 



-1^0, 



(40) 



for g — > oo. We might validate this condition for a 
2DEG with the potential of the inter-electron interaction 
Vee{q) oc 1/9 by using the fact that the polarizability 

I 



function x('z) ~^ for g — >■ 00^ Actually, even for finite 
x(q), as obtained in the long- wavelength limit— or for 
quantizing magnetic fields p^ when the left-hand side of 
Eq. (|40p is finite, the integral over the path p2 is con- 
vergent and falls off with distance r without oscillations. 
The integral over the path p^ also has this property. 

As a consequence, the oscillatory behavior of the right- 
hand side of Eq. (P^ is entirely defined by the poles of the 
function 1/e (g, il) located inside the integration contour. 
To find them, we can use Eq. (|T5|) . The poles coincide 
with those of the function l/[fi — LOp{q) + i/r] because of 
the adoption fi > and ujp{q) > 0. Again, calculating 
the terms with n = — 1 in Eqs. ([TO)) and (pHl) . we deal 
with the poles of the function l/[— ^ -I- ujp{qx,0) + i/r] 
only. 



Repeating the reasoning of Sec. IIVI we suppose that 



the poles q, = q'^ + iq'- (j = 1, 2, 
small imaginary parts: \q','\ ^ 



.) are simple and have 
'■|. For the real parts 



q'j we have LOp{q'j) — 51, while for the imaginary parts 
q'J = l/uj'p{q'j)T, where uj'p{q) = dujp{q)/dq. The poles 
are inside the integration contour if lo' (q'A > 0. Keeping 
only the oscillatory contributions from the poles, we have 
for cos (6 — 9) > X 



$+ = $(r,0,, 



nn 



s{4,-e)>\ 



E 



j{<>0) 



q'jVo{q',)Ai{q'^ COS cb,q'^ sin (j)) ,(^,+,,.),,„,(^_,) 



(41) 



Analogously, for tt/2 + 9 + X < (j) < 3tt/2 + 6 — X, the original integration path in Eq. (jJH) should be supplemented 
by the quarter-circle path l^l = const -^ -|-cx) with — 7r/2 < argg < and by the path over the imaginary axis Imq < 
to form a contour. The plasmon poles are inside the contour of integration if uj' (q'A < 0. The result is 



$_ = $(r, ( 



s{<p-g)<-x 



iTrft y 



q'jVoiq'j)Ai{q'^ cos <j,,qr sin <P) ,(,^.+,,y),eos(0-e) 



i('^;<o) 
Returning to Eq. (j38|) . we obtain the following expression: 



u^'piq-) 



r(r,6i)- / <P+d(j)+ I $_, 

-7r/2+e+A tt/2+B+\ 



O $ 



Ai 



7r/2+6l 



37r/2+e 



^+d4) 



$_ d0 + O $ 



Ai 



-7r/2^ 



it/2+0 



(42) 



(43) 



r 



so that we have an asymptotic representation of the func- 
tion $ (r, 9, (j)) for the whole domain of the variable ip. 

The integrals in Eq. (|43|) can be evaluated with the help function H,„ (y) -^ 
of the approximation 



which are obtained from integral representations of the 
Bessel function of the first kind Jm (y) and the Struve 



Ai (qcoscj), qsint 
and the integrals 



-iq{X cos (f> + Y sine/)) , (44) 



sin (y sin x) 
cos (y sin a;) 



dx 



Joiy) 



(46) 



smx 



sin (y sin a;) 1 ^^ ^ ^ f Ji{y) \ 
cos(ysinx)J \H_i (y)/ 



supplemented by derivative identities 



(45) 



d-y^' 



Jmiy) 
H™(y) 



= y 



Jm-i (y) 
Hm-i (y) 



(47) 
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and the property J-m{y) = (— l)'"<^m(y) (the latter is 
vaHd for integer m). To evaluate the potential of the 
frozen wave, given by Eq. (|20|) . we also use the following 
result deduced from Eqs. (|^ and (|T7|) : 



cos (y sm x) 



dx = IT 



H-i(y) 

y 
Jiiv) 



~ H-2 (y) 

^2 (y) . 



(48) 
Having established the integration procedure, we 
merely present the result for the right-hand sides of 
Eqs. P^ and (1^ . We suppose that Q > is a single 
solution to the equation ui^iff) = fi so that co'piQ) > 
(see Fig. [T]). For circular polarizations of the MWs, 
£y =^ H£x, the principal wave is a circular plasmon coro- 

I 



fating with the electric field vector of the radiation [coun- 
terrotating if Wp(Q) < 0]: 



V^^\r,9,t) ~ -^(e±'''-'"*+*'^i^i(r) +C.C.), (49) 

where the undamped amplitude of the wave is 

27:^Vo{Q)QMX\ 



Vd 



<{Q) 



(50) 



c.c. stands for the complex conjugate of the preceding 
term in parentheses of Eq. (^^1) . and, as a reminder, %j) = 
argX. The complex radial function Fi{r) is 



i^i (r) = iJi 



^'p{Q)T 



+ H_i Q 



^'piQ)r 



Oi^ 



(51) 



I 

Let us inspect the asymptotics of the Bessel and Struve for the function Fi (r) : 
functions'^^ 



Jrn (0 






(52) 



Fi(r) 



TvQr 



iQr-r/uj' (Q)T-iT^/i 



(54) 



H™(e) = \/^sin(e 
V 7rf V 



~~ 4 



as Qr — > oo. Note that we have used Q ^ 1/uj' {Q)t 
under the square root in the expression above. 

For the frozen plasmon induced by MWs with circular 
O (C™'^) I (53) polarization we have the following isotropic potential: 



VdQ\X\ 



Fo{r)+c.c.y 



valid for ^ — > oo [below we use Hm (0 only for m < 0, so ^ v) 

that the form of Eq. ([55)1 is adequate for our goals] . Using 

Eqs. ([52]) and ([55)1 in Eq. (ISTt . we obtain the asymptotics where the complex radial function Fo{r) is 

I 



(55) 



Fo (r) =iJo[Qr + — 



^piQ)^ 



- Ho Q 



o'piQ)' 



o 



Ao 



(56) 



valid for the far zone r ^ Aq (see Appendix [C|). The 
function Fq (r) has the asymptotics, similar to -Fi (r) 



Fo{r) 



^iQr-r/Lj'^{Q)T-iTT/4 



(57) 



for Qr — > oo. When compared to the principal plasmon, 
the amplitude of the frozen wake has an extra small pa- 
rameter Q\X\ as a factor, similar to the wire defect case 
[seeEq. dag. 

Let us obtain the potentials V''^'^{r,d,t) and V^°'^{r,6) 
for the case of a linear polarization of MWs, for example, 
for X ^ and Y = (without magnetic field this is 



I 

equivalent to £y =0). We obtain the principal plasmon 
having the character of a standing wave: 

V'^^\r,0,t) - -^(e-*"*+*^i^i(r) +C.C.) cos 9, (58) 

with the maximum amplitude along the x axis, 9 = 0. 
For the frozen plasmon we have the following result: 

F(") (r, 9)^- ^£^M (^Fo (r) sin^ 9 + F2 (r) cos 20 + c.c.) , 

(59) 
where the complex radial function F2{r) is 



F2ir) 



Ji(Qr 



<(Q)r 



iJ2[Q 



t^'piQ)r 



H_i gr 



<iQ)r 



Qr 



H 



Qr 



^'p{Q)r 



o\^ 



, (60) 



which has the asymptotics F2 (r) ^ Fq (r), so that m the 
hmit Qr — > 00 we have the fohowing result for the frozen 
wave: 

V(°\r,e) ^ -^^^M^(^o(r) +C.C.) cos^ 0, (61) 

where for the function Fo{r) one should use Eq. ([57|. 



I 

set that again has the fundamental lattice period 2tt/Q. 
Due to the fact that even without damping the plasma 
wave from a dot defect decays with distance as (x Xj \JQr^ 
the interference pattern has relatively low contrast. To 
demonstrate this, in Figs. U and [S] we plot the function 



N-\ 



Vf^a;, y) = ^ 1/(0) (x, y - 2^nA/Q) 



(62) 




FIG. 4. Potential (in units of PdQ|X|) of the net frozen 
plasmon from A'' = 10 dots, which form a truncated periodic 
set occupying the space 2: = 0, < j/ < (A'^ — 1) (in units 
of 2-k/Q). Constructive interference. Circular polarization. 
Damping is neglected. 




FIG. 5. Potential (in units of 'PdQ\X\) of the net frozen plas- 
mon from N = 10 dots, which form a truncated periodic set 
occupying the space x = Q, Q < y < {N — l)/2 (in units 
of 2tv/Q). Destructive interference. Circular polarization. 
Damping is neglected. 

Finally, to amplify the plasmon wave with a set of dot 
defects, we should arrange them in a truncated periodic 



which represents the frozen plasmon wave produced by a 
linear chain of N dots separated by the distance 27rA/(5, 
for two values A = 1 and A = 1/2 for circular polariza- 
tion of the radiation. 



VII. ROLE OF RETARDATION AND 
POLARIZATION OF DEFECTS 

We have neglected the retardation of the electromag- 
netic interaction, focusing on plasma excitations with the 
wave numbers Q ^ k, where k is the wave number of the 
MW radiation. This has provided us with the approxi- 
mation that the spectra of plasmons in the rest and the 
moving reference frames coincide with each other (be- 
cause Poisson's equation is invariant with respect to a 
uniform coordinate transformation). The plasmons thus 
obtained do not experience the radiative decay?^ the ef- 
fect being recovered by the inclusion of the retardatiouj^^ 
In a clean 2DEG, this makes the excitation of plasmons 
with larger wave numbers even more efficient because 
they are less prone to decay. 

The retardation can be implemented in the developed 
method by following the line of the theory of interaction 
of relativistic particles with matter, ^^^ although it is com- 
plicated by the necessity to allow also for the spatial dis- 
persion of the conductivity tensor. ^^ Note that Maxwell's 
equations change in the moving reference frame due to 
the non-Lorentz-invariant transformation realized by the 
operator S. The excited plasmons can still be obtained 
analytically for wave numbers Q <C 1/R. On the other 
hand, the existing approaches j^^— which are devoid of 
the relativistic limitations, would require extensive nu- 
merical simulations if they were modified to take into 
account collective excitations with large wave numbers. 
Such are in particular plasmons with negative group 
velocity, which produce backward propagating plasmon 
waves, readily obtained within the wake- field method as 
long as their spectrum LOp{q) and the relaxation time t 
are known. 

A more delicate issue is the radiation-induced polar- 
ization of the defects, leading to the scattering of the 
incoming MWs. We have neglected it, regarding the de- 
fects as elementary particles lacking an internal structure. 
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Strictly speaking, this approximation is valid only for de- 
fects that are represented by quantum mechanical objects 
with the excited states lying far above the energy of the 
radiation, and in weak radiation fields so that the dipole 
moment of the defects multiplied by the electric field is 
small as compared to the energy of the excited states. 
This is not a typical situation. On the contrary, the 
defects are usually classical objects: ripplonSj^ lattices 
of metal scatterers and gratings^i^ or surface acoustic 
wavesi^ 

To make allowance for the general situation, wc should 
first determine the potential V^^'^'^'* (r, t) of the scattered 
wave and carry on with its sum with the potential of the 



bare defect forming the total unscreened potential Vq 



(tot). 



Vr 



(tot) 



(scat) 



ir,t)^Voir) + Vr^"ir,t) 



(63) 



The scattered wave contains frequency harmonics fls, dis- 
crete for brevity's sake, which may include those not 
present in the external field E(t): 



yo(^^^*)(r,t) = ^y„(^^^')(r,a)e^ 



ifl^t 



(64) 



o. 



It suffices to mention a two-level system whose dipole 
moment oscillates not only at the frequencies fig = ±f^ 
of the external field but also at the frequencies fig = 
±f2 ± Or, where Or is the generalized Rabi frequencyi^ 
In the moving reference frame we naturally have 



Vo^'°'\r,t) = S^V^'°'\r)S 



= Voir-ro{t)) + vt^'\r-roit),t) 



(65) 



This potential is then screened by the 2DEG and pro- 
cessed in the straightforward way developed above. 

For the principal plasma wave, in the linear regime, 
we can neglect the weak MW-induced polarization of the 
2DEG in the scattered wave potential. This is equivalent 
to the following approximation: 



Vr 



(tot) 



{r,t)^Vo{r-ro{t)) + V^''^'^{r,t) 



(66) 



In this approximation, the excited plasma wave consists 
of two contributions: the one originating due to the 
MW-induced polarization of the electron gas and scat- 
tering of these polarization waves at bare defects, con- 
sidered in detail in this paper, and the contribution that 
is the MW-induced polarization potential of the defects 
screened by the MW-unpolarized electron gas. The lat- 
ter is sometimes treated as the only source of the excited 
plasmonS)^ which is reasonable if the defects are not 
charged and located far enough from the electron gas so 
that they weakly influence in particular the dc conduc- 
tivity of the electron gas. These two excitation mecha- 
nisms cannot be separated in the local permittivity- or 
conductivity-based approaches J^"— 

To complete the analysis, let us compare the strengths 
of these two contributions for a defect given by a metal 



sphere with Coulomb charge Z, radius Rq, and center 
with the coordinates (x — 0, y — 0, z — d), isolated from 
a 2DEG at z = 0, so that d > Rq. The system is placed 
in a medium with permittivity k. The length d is small 
as compared to the wavelength of the MW radiation in 
the medium 2ttc/Q^/k, where c is the speed of light in 
vacuum. Under these conditions, the 2D potential of the 
defect is 

Vo{q) = -^c-^'. (67) 

ZTTKq 

The 3D polarization field in the conditions of small dis- 



tances ii 



,41 



Vr 



(scat 3D) 



(r,z,i) 



i?3r • E(t) 



(r2 + {z- d)2)3/2 



(68) 



[we have previously defined r and E(i) as 2D vectors], 
which corresponds to the polarization of a sphere in a 
static homogeneous electric fieldi^ The spatial Fourier 
transform of the 2D scattered wave potential is thus 



Vr 



(--t)(q,t)^!Me-9'^q.E(i). 
27rq 



(69) 



This potential screened by the 2DEG and reduced by 
its bare value has the following form in the space-time 
coordinates: 






-1 



(fe4")-t-g,4"))e^'i'-™"*d2^. 



(70) 



£j and£J"^^ ^£*. 



where £" ■ ' 

The structure of Eq. ^ is identical to that of Eq. ^ 
in which the potential is expressed by Eq. (p7|) . More- 
over, without magnetic field, X (x £x and y ex £y, so that 
the two excitation mechanisms produce additive contri- 
butions, with the relative sign being determined by the 
sign of the charge Z [see Eqs. ([6]) and (|44| ]. In finite mag- 
netic fields these waves interfere. Let us finally compare 
their relative strength, at Wc = for definiteness. We 
introduce the electric potential of the sphere II = Z/kRq 
and the wavelength of the MW radiation in vacuum, 
I — 2ttc/Q, and analyze the following dimensionless pa- 
rameter: 



71 



en 



/ 



471^171*0 \Ro 



(71) 



If 7?. ::!> 1, the polarization of the defect can be neglected 
even though it is a classical object. To check whether 
this is reasonably achievable in a GaAs/AlGaAs-based 
2DEG, we use the following benchmark data: decent 



ell = 1 eV while An^ 



1.4 MeV. So we require 



that Ro <$: I X lO"'^, and if Rq 



1 fj,Ta, the "ele- 



mentary particles" mechanism dominates for frequencies 
il/2TT ^ 0.3 THz. This figure is somewhat below the 
range of interest for actual applications,* yet there is a 
definite reserve in increase of the potential II and decrease 
of the size of the defects Rq . 
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VIII. CONCLUSIONS 

We have developed an efRcient method to find plas- 
mons excited by MWs with a homogeneous field in a 
2DEG with parabolic spectrum having defects, treated 
as weak scatterers. The method uses the concept of the 
wake created by a charged particle moving in plasma. 
The expressions that we have derived for excited plas- 
mons explicitly address the efficiency of the photon-to- 
plasmon conversion and the type of the resulting excita- 
tions. 

We have shown that the photon-to-plasmon conversion 
is more efficient for plasmons with larger wave numbers, 
even though they might not be well sustained in the sys- 
tem due to damping. Intense persistent plasma waves 
are created at optimal wave numbers that depend on the 
spectrum of plasmons and the distance that the waves 
are required to travel without being substantially atten- 
uated. 

Concerning the type of excited plasmons, they can 
have traveling or standing wave components, depending 
on the geometry of the defects and the polarization of 
the radiation. Thus, a single straight wire defect can 
produce traverse plasmons traveling away from the wire 
if their group velocity is positive [see Eq. (1^51) ] and toward 
the wire if it is negative [see Appendix [B] . A single dot 
defect can excite circular plasmons [see Eq. (|49|] when 
the MW radiation has circular polarization, or standing 
plasmon waves for linearly polarized MWs [see Eq. ([55]) ] . 

I 



The zeroth harmonic of the excited plasmons, which is a 
frozen charge density wave, has also been obtained [see 
Eq. (|29|) for a wire defect and Eqs. (f55|) and (f59|) for a 
dot defect]. Defects arranged in truncated periodic sets 
can create amplified interference patterns, which can be 
found easily [see Eqs. (I55|) and (|36p and Fig. [5] for wire 
defects, and Fig. |3] for dot defects]. 
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Appendix A: Definition of w^ake far zone for wire 
defect 



If we specify the boundaries of the wire as \x\ = a, for 
the Fourier transform Vw{qx) of its potential Vwix) we 
can use the following expression: 



zvr 



(Al) 



Let us rewrite Eq. ([2T|) explicitly using Eqs. (|TOl) and 
dAi]) for Al {qx,0) and Vn,{qx)-- 



I{x) 



n 



{2^Y 



7r/n a -f-oo 

dte'^^ I dx'V^{x') I dqx 



1 



e(9x,0,f7) 



JqJx-x' -xa{t)\ 



(A2) 



The expression in square brackets in Eq. (jASj) is posi- 
tive if a; > a -I- 2|A| and negative if a; < —a — 2\X\ for 
the entire domain spanned by the integration variables: 
—a < x' < a and —n/fl < t < ir/n. Thus, we define 
the far zone, jxj > a -I- 2|A|, as the region where struc- 
tural details of the functions Ai {qx,0) and Vw{qx) do 
not influence the shape of the wake. The value obtained 
for the position of the far zone boundary is applicable 
to the principal plasmon given by Eq. (jl9|) . The right- 
hand side of Eq. (pn|) has an extra factor An under the 
integral, which results in an extension of the far zone to 
|a;| > a -(- 4|A| for the frozen wave V^^\x). 



Appendix B: Excitation of plasmons with arbitrary 
spectrum 

Let us briefly discuss different situations that may 
arise while considering collective excitations with arbi- 



I 

trary spectra. As an example, in Fig. [5] we present the 
spectrum of the principal 2D magnetoplasmon mode in 
a weak magnetic field (several Landau levels are filled), 
obtained in the random phase approximation)^ Nega- 
tive qx are not shown, but inferred. For qx > with 
qy = 0, the spectrum and the radiation harmonic uj = fl 
have three common points. The case corresponding to 
the point A is discussed in the main text. The point B 
(at q^ = q[) is treated analogously to the point A, but 
here the radiation excites plasmons with negative group 
velocity uj'p{q'i,0) < at positive qx — q'l with a com- 
plementary uj' {—q[,0) > at negative qx = —q'l (not 
shown). This generates a component of the principal 
plasmon wave, analogous to the one given by Eq. ([26t . 
traveling toward the wire defect with the phase velocity 
Vp (x) — — sgn(x) i^/q[ (neglecting damping). 

The point C is a tangency point, for which uj' [qi^^ 0) — 
0. The approach used in the main text is not valid now. 
Actually, considering the model integral of Eq. (HI]), we 
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0), a.u 




q , a.u. 



FIG. 6. Spectrum of the principal 2D magnetoplasmon mode 
in weak magnetic fields as a function of q^^ at Qy = (full line) 
and the radiation harmonic lj — il (dashed line). Negative 
qx are not shown. Contributions of the plasmons with the 
wave numbers q^ — q'l (the intersection point B) and qx = q'2 
(point C) are discussed in Appendix B. Interaction of photons 
with plasmons at smaller wave numbers (point A) is detailed 
in the main text. 



have Sq^ « 21/0;^ (gj, 0)t, which gives two complex values 
of 5q. The poles of the integrand in Eq. ()21|) remain sim- 
ple and located at Qx — q2 + ^1- Only one of them falls 
inside the integration contour for positive x, the other 
contributing for negative x. This gives two extra plasma 
waves (totaling four with points A and B) propagating 
in opposite directions: the one traveling away from the 
wire has the wave number (5_|_, while the wave traveling 
toward the wire has the wave number (3„, where 



Q±=g^±sgn(^;'(g^,0)) 



1 



K'(9^,0)|r- 



(B3) 



Both waves are characterized by the exponential decre- 
ment exp (-|a;|/^|a;;;(<7^,0)|r| 



must retreat to find the roots of the equation 
Wp(fc,0) = Vl + i/T. 



(Bl) 



For weak damping, we can use the next term in the Taylor 
series expansion around the point gj- 



Appendix C: Wake far zone for dot defect 

If the dot defect occupies the area ttD^, where _D > 
is the nominal radius of the defect, for the Fourier trans- 
form Vo(q) of its potential Vo(r) = Vo(a;, y) we can use 
the following expression: 



c.p(fe,0)«c.p(g^,0)+c.;(g^,0)5g + c.;'(g^,0)^, (B2) V^{q) ^ // l^li^ e-''(-'-^'^+^'™^) da^W, (CI) 



where 6q ^ q^ — q'2 and ujp{q2,0) — Q. The expansion 
given by Eq. ()B2p can be used when the third term in 
the right-hand side is not small as compared to the sec- 
ond one, in particular, at the tangency point C where 
ui'p{q'2,0) = 0. In this case from Eqs. ((BT|) and ((B2|) we 



r'<D 



{2^r 



where we have used polar coordinates for the reciprocal 
space, as in Eq. ([39]). and have kept the superfluous de- 
pendence on the angular coordinate in the integrand. 
Let us rearrange Eq. ([55]) explicitly using Eqs. (flU)) and 
(IC1|) for Ai (q cos 0, g sin 0) and Vo{q): 



$ 



^ 



(2^r 



tt/O 



dte 



int 



-TT/n 



+00 



dx'dy'Vo{x',y') 

r'<D 



dqq 



eiq,n) 



^igp 



(C2) 



where the scalar quantity p, which is a function of the 
variables {x',y',t) and the parameters {r,9,(j)), is 

p = r cos {(f) — 9) — [x' + xq (t)] cos — [y' + j/o (i)] sin (f>. 

(C3) 
Our goal is to find conditions when this operator is posi- 
tive and negative definite. Such conditions must depend 
neither on the variable t, which characterizes the MW- 
related coordinates xq (t) and yo (t) j nor on the variables 
x' and y', which specify the structure of the defect. Then, 
as in the case of a wire defect, the shape of the wake will 
not depend on structural details of the functions Ai (q) 
and Vo{q). 

In contrast to the wire defect case, the boundary of the 
far zone can not be rigorously defined. Let us introduce 



r 



the following quantity; 



A, = J{D + 2\X\f + iD + 2\Y\f. (C4) 



Then, for r ^ Ai we have 

f/5>0, -f + A< 
1/5<0, f + A<(A 



<f-A, 



< 



Stt 



A, 



(C5) 



where A « Ai/r — > as r — > 00. 

The expression that we obtain in the main text for 
the plasma wave created by a dot defect is valid in the 
asymptotic limit r/Ai — )■ 00. 

For the frozen plasmon, the right-hand side of Eq. ([20)1 
has an extra factor A„ under the integral, which makes 
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us introduce 



Ao = VP + 4|X|f + (D + 4|y|)^ (C6) 



and deal with V^^'^ (r) in the asymptotic hmit r/Ag -^ oo. 
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